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ABSTRACT 


Modeling and simulation have proven invaluable to the design of various systems 
and planning of many experiments since the free-electron laser was invented in the 
1970s. This thesis illustrates methods used for these tasks and for the validation and 
development of further theory. FEL systems and their components are described, briefly 
discussing how they fit into design considerations. Basic theory of FEL operation is 
reviewed, including the resonance condition, FEL wave and pendulum equations, and the 
concepts of undulator taper, extraction and beam quality. Two computer programs used 
for FEL simulation (FEL3D/FEL4D, developed at the Naval Postgraduate School, and 
Ginger, developed at Lawrence Berkeley National Laboratory) are discussed and 
compared, then applied to specific FEL cases. The cases studied serve to evaluate key 
performance characteristics of FELs. The final case study directly compares output from 


the two programs. 


In addition to this type of physics simulation, computer tools are also vital for 
component design. Relevant methods are illustrated using the example of an electrostatic 
cathode test cell and drawings made using the Los Alamos National Laboratories design 
tool, Poisson Superfish. After reviewing the principles involved and system constraints, 
a new test cell is designed for use with future NPS experiments. The effects of anode- 
cathode spacing and anode aperture size are examined for two different anode 


configurations and recommendations are made for future design iterations. 
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I. INTRODUCTION 


Free-electron lasers (FELs) are, like any other lasers, sources of coherent 
electromagnetic radiation (light), which is amplified by stimulated emission from 
electrons. Free-electron lasers differ from other types, however, in how the light is 
generated. Rather than relying on electrons bound within matter, FELs use a high-energy 
beam of electrons in vacuum. They are distinguished from solid-state and gas lasers by 
the flexibility of their designs: Adjustment of a few system parameters completely alters 
their output wavelength, whereas conventional lasers are restricted by the quantum 
properties of their lasing media. Some of the relevant FEL parameters will be set during 
system construction; but others, such as the electron beam energy, can be altered 


continuously, allowing wide flexibility in application. 


Years of study with FELs have also shown the systems to be powerful, reliable 
(98% or more uptime on many systems) and efficient (10%—20% wall-plug efficiency 
according to some estimates). On the other hand, FEL systems have historically been 
large and expensive; a smaller footprint and the ability to be supported on mobile 


platforms are active areas of research for some applications. 


Modeling and simulation have proven invaluable to the design of various systems 
and planning of many experiments since the free-electron laser was invented in the 
1970s. This thesis illustrates methods used for these tasks and for the validation and 


development of further theory. 


Chapter II describes FEL systems and their components, briefly discussing how 
they fit into design considerations. Chapter III reviews the basic theory of FEL 
operation, including the resonance condition, FEL wave and pendulum equations, and the 
concepts of undulator taper, extraction and beam quality. Chapter IV discusses and 
compares two computer programs used for FEL simulation (FEL3D/FEL4D, developed 
at the Naval Postgraduate School, and Ginger, developed at Lawrence Berkeley National 


Laboratory) and applies them to study specific FEL cases. The cases studied serve to 


familiarize readers with the simulations’ output and highlight key performance 
characteristics of FELs. The final case study directly compares output from the two 


programs. 


Chapter V delves into the specifics of physical design, using the example of an 
electrostatic cathode test cell and drawings made using the Los Alamos National 
Laboratories design tool, Poisson Superfish. The principles and design constraints for a 
test cell to be used with future NPS experiments are examined. The effects of anode- 
cathode spacing and anode aperture size are studied for two different anode 


configurations and recommendations are made for future design iterations. 


Il. KFREE-ELECTRON LASER SYSTEM DESCRIPTION 


Free-electron lasers can have two basic configurations: oscillator or amplifier. 
Most components are common to both designs, as shown in Figure 1. The difference is 
how the light is extracted from the electron beam—amplifiers amplify an existing optical 
field over one pass while oscillators use optical resonators similar to conventional lasers. 
The following sections discuss the purpose and general properties of the system 


components. Figure | depicts a particular FEL system with a few illustrative properties: 


1. Separate beam paths for running in oscillator or amplifier mode. This 
shows the difference between the two, but the dual beamline configuration can be found 


in real systems, such as the experimental setup at Thomas Jefferson National Laboratory. 


2. A superconducting radio-frequency linear accelerator, and the helium 
refrigeration plant associated with it. Normal-conducting systems are also used widely. 


Electrostatic (DC) and storage ring accelerators have also been used for FELs. 


3. Beam recirculation for energy recovery (see “Beam Disposal” below). 








(ie aM 1+ Pi pa 
Undulators yt 4 4 
FEL Amplifier 


Figure 1. FEL System Schematic 





A. INJECTOR 


The injector generates the beam of free electrons from a cathode. Cathodes come 
in several varieties and generally produce electrons either via thermionic emission, the 
photoelectric effect or field emission. Thermionic cathodes emit electrons when 
electrically heated to a point where the thermal energy in the metal exceeds its work 
function. Photocathodes emit electrons due to the photoelectric effect. When they are 
illuminated by a drive laser, as long as the energy per photon exceeds the work function, 
electrons are liberated. Field-emission cathodes are not necessarily metal, and use 
geometry to create points where high electric fields exist once an external potential is 
applied. For example, diamond field-emitter array (DFEA) cathodes are wafers of 
diamond film grown on a suitable substrate such that the diamond forms an array of 
microscopic pyramids; the field is enhanced enough at the pyramids’ tips that electrons 


escape by Folwer-Nordheim tunneling. 


Regardless of cathode type, injectors will apply some external field to the emitted 
electrons, accelerating them to modest energies before they enter the accelerator 
component. There will also be some amount of instrumentation and beam transport 
apparatus (e.g., focusing magnets) to prepare the beam for entry to the accelerator and to 


perform initial characterization of the electron beam. 


Injector design can be important to overall system performance in that it must 
provide electrons to the accelerator with optimal bunch length and repetition frequency. 
Furthermore, the beam quality (i.e., the spatial or angular spread of the electrons or the 
relative disparity of their kinetic energies; see discussion in Chapter III) provided by the 


injector is the highest that will be present in the entire system. 


B. ACCELERATOR 


The electron beam is next passed through an accelerator, which raises its energy 
to the desired level. Direct-current (DC) acceleration is possible, but high-energy 
applications usually use radio-frequency (RF) fields. Electron accelerators generally 


consist of a series of evacuated metal cavities to which external RF power is applied in 


order to generate strong electromagnetic fields inside. The alternating fields can be timed 
so as to provide net positive acceleration to the sequence of electron pulses provided by 


the injector. 


C. UNDULATOR 


The undulator consists of a series of alternating transverse magnetic fields, which 
produce transverse accelerations of the electron beam. This causes the electrons to 
radiate, producing spontaneous emission like population inversion does in a conventional 
laser. The transverse motion of the electrons also allows for transfer of energy between 
the electrons and the optical field, stimulating further emission, creating a coherent light 
source. The relativistic speeds of the electrons ensure this optical field takes the form of 


a narrow coaxial beam. 


Undulators can be linear or helical in their polarization. Linear undulators (see 
Figure 2) consist of a series of dipole magnets, alternating in orientation, which create a 


magnetic a field that is periodic in one transverse plane. 





Figure 2. Linear Undulator, showing (1) Magnets, (2) Electron Beam and (3) 
Synchrotron Radiation (not discussed) [After 1] 


Helical undulators, as shown in Figure 3, usually consist of permanent magnets or 
electromagnets coiled around the beam path in a helix, and provide a magnetic field that 
is periodic in both transverse directions (see Chapter III). There are also planar-geometry 


undulators that can provide a helical magnetic field, such as the APPLE (Advanced 


Planar Polarized Light Emitter) system [2]. In the APPLE, two linear undulators are 
operated with their planes of polarization at an angle to one another, thus producing net 
polarization that can be helical (elliptical or circular) or linear in a specified direction, 


depending to the angle used. 





Beam Pipe: 1D =5.3mm 
Low-Carbon Stee/ Pole OD=6.3mm 


Figure 3. Helical Undulator [From 3] 


D. OPTICAL COMPONENTS 


In FEL oscillator designs, the undulator is contained in a resonator cavity, with 
one fully reflective and one partially transmissive mirror. As the light repeatedly transits 
this resonator, it serves to stimulate further emission from new electrons. On each pass, a 
portion of the light is outcoupled via the partially transmissive mirror. In most FEL 
amplifier designs, the initial optical field is produced by a seed laser that is amplified 
while passing only once through the undulator. There are also self-amplifying 
stimulated-emission (SASE) designs that use the spontaneously-emitted light from the 
first portion of the undulator as a seed. Amplifier designs, in general, require longer 


undulators than oscillators, with SASEs requiring significantly more length. 


E. BEAM TRANSPORT AND DIAGNOSTICS 


The physical layout of the FEL system will determine exactly which components, 


such as bending magnets and focusing magnets, are required to properly transport the 


beam. Most systems will also have several types of diagnostics at various places along 
the beamline, both for application-specific purposes and to track the results of beam 


transport. 


F. BEAM DISPOSAL 


After light is generated in the undulator, one must dispose of the spent electrons. 
Ultimately, the beam’s energy can be dissipated into something as simple as a metal 
target as long as appropriate cooling and radiation shielding are provided. This is labeled 
as “Beam Dump” in Figure 1. Since light generation only depletes a small fraction of the 
electrons beam’s energy, some designs recirculate the beam through the accelerator; if re- 
injected with the proper phase, the electrons will give up significant energy to the RF 
fields, boosting system efficiency and making beam disposal much less hazardous. This 
technique is referred to as “energy recovery,” and linear accelerators that use it (not all of 


which are FELs) are known as “energy-recovery linacs” or ERLs. 
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Ul. BASIC FREE ELECTRON LASER THEORY 


The centerpiece of FEL physics is the interaction between electrons and light that 
takes place within the undulator. The physics of particle acceleration and laser beam 
transport are vital to the operation of the greater system; but the undulator is where the 


interaction characteristics of free-electron lasers are determined. 


It is helpful to begin a mathematical description of these interactions by defining a 
few useful coordinates: 


ct 


¢ =(k+k,)z—ot, and {3.1} 


pe eee ye ai 
dt 


The normalized time 7 serves as a temporal variable along the longitudinal spatial 
coordinate; t = 0 corresponds to the beginning of the undulator and 7 =1 to its end. The 


electron phase ¢ measures an electron’s relative position in a portion of the electron 


beam the length of one optical wavelength 4, and the electron phase velocity v 
measures its rate of change. Here z is the electron’s position along the undulator; 
@=kcis the frequency of the optical field; L is the length of the undulator; 


k=2nr/Aand k,=2a/A, are the wave numbers corresponding to the optical 
wavelength 2 and undulator period A, respectively; and £, =v_,/c is the normalized 


longitudinal electron velocity. We will see that these definitions play an important role in 


the free-electron laser’s interaction. 


A. RESONANCE CONDITION AND PENDULUM EQUATION 


The electrons’ motion is determined by the Lorentz force equations, namely 


408) __© (84 8x8), and 


dt mc 
3.2 
dy__&@ 3 ve 
dt mc 


For a helical undulator (see Chapter II), the magnetic field can be expressed as 


B,,, = B(cos[kyz], sin[k,z], 0), {3.3} 
while the optical fields will take the form 

E = E(cosy,—siny,0) 

B,, = E(siny, cosy,0), where {3.4} 


y=kz-att+@. 
In the preceding equations, EF is the optical electric field strength, B is the undulator’s 
magnetic field strength, @ is the field phase, and e and m are the charge magnitude and 
mass of an electron. In Equation {3.2}, the pertinent magnetic field is the sum 


B=B_, +B. 


Substituting {3.3} and {3.4} into {3.2} and assuming that f£, ~1 gives the 


transverse motion as 


a: OS 


B = (cos[k,z],sin[k,z],0), where 
Y 
13.5} 
K = CB nshy . 
2amc* 


The “undulator parameter” K is a dimensionless measure of field strength and is of order 
unity. Note that this derivation assumes ideal injection into a helical undulator such that, 
when solving for position as a function of time, the constants of integration are all zero. 


In the helical case, B,. = B; but a similar analysis holds for linear undulator geometry 
with B,. =B/ 2. 


The evolution of electron energy can be analyzed through changes in y and in v. 
In particular, from the second equation in {3.2} and Equation {3.5}: 


EE ceed {3.6} 
dt ymc 


The electron’s velocity is related to its energy and hence its Lorentz factor y, specifically, 


y” =1- fi -£. {3.7} 


10 


Equation {3.5} shows that the magnitude of the transverse velocity is 


merely K /y ; thus, for vy >> 1, Equation {3.7} reduces to 


2 2 
TO ell pee 3.8} 
y 2y 


As an aside, we consider next the condition under which optimum energy transfer 





occurs between the electron beam and optical field—treferred to as “resonance.” In this 
state, each electron emits (or amplifies) one wavelength’s worth of light (42) for each 
undulator period (distance /,) through which it passes. If the interaction between the 
electron and the optical field is thought of as a race down the undulator, then the light 


leads in velocity by c(1— £,) and wins by a distance of 2 in one undulator period. The 


time passed is 








A 
At=—. {3.9} 
cB. 
The winning distance can be expressed as the velocity difference multiplied by the race 
time: 
i 
A=c(l— B,)At=A, B, {3.10} 
B. 
From Equation {3.8}, we see that the condition for resonance is 
4.(4-1}, or, for y > 1, 
0 : {3.11} 
REA oe 
27 


which relates the FEL’s output wavelength to fundamental properties of the electron 
beam and undulator. This important relation highlights the tunable and designable nature 
of FELs. Adjustment of any or all of the system parameters that impact K allows 
designers to build systems that emit at any of a wide range of wavelengths, or even to 


adjust K on an existing system. 


Using the “resonance condition,” Equation {3.11}, in the definition of phase 


velocity, and recognizing that k >> k, (i.e., vy >> 1) shows that 
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{3.12} 





From the definition of undulator period (ZL = N /,, where N is the number of periods), 


we find 


Azan 3.13} 
y 


Differentiating v with respect to t and using the result of Equation {3.6} gives: 





2 
ae = = |a|cos(¢ + ~), where 
dt dt {3.14} 
_ 4aNeKLE 
lal = Pix i20 
y°me 


This is known as the “FEL pendulum equation” for its similarities to the equation of 
motion for a simple pendulum. Using this equation, the evolution of FEL electrons is 


often described with phase space (¢ vs. v) graphs. 


Parameter a is a dimensionless measure of the laser field amplitude EF. Values of 
a<q are considered the “weak field” regime, while a>>z is considered “strong 
fields.” The FEL theory as presented is valid for both regimes, as well as intermediate 


values. 


B. FEL WAVE EQUATION 


The FEL’s optical field will evolve according to Maxwell’s wave equation: 


fv 253 }Ae. = “2 Fea. {3.15} 
Cc Cc 


In this equation, A is the optical vector potential and J is the transverse current density 
due to the undulations of the electrons. The vector potential is related to the optical fields 


by: 


B=---=, B-=vVxA. {3.16} 
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Assuming a complex optical electric field in the form 


E (x,t) =|E(%, |e?" e*, {3.17} 





where a = kz — at represents the carrier wave, the vector potential becomes 


A(z, 0) = ie", {3.183 


where 7 is the polarization vector. Given that the light spectrum produced is not only 
finite but narrow (i.e., a laser), it is safe to assume that the amplitude and phase vary only 
slowly along the scale of an optical wavelength (4 = 2z/k ) and the scale of the optical 
frequency (@=kc ). In this case, one can assert 


: OO ea Fl cle 
Ot 0 


o\z| 
—_ <« @|E at 
Ot Zz 








= <ko. 3.19} 


ae 
Substituting Equation {3.18} into Equation {3.15} and using conditions {3.19} to 


neglect terms with two derivatives, the relation simplifies to 


i>. 10 7 Cd ae ere 
[vi +2i( 2422) a3. --SA (7-0 )e {3.20} 


where the transverse Laplacian symbol (V7) denotes differentiation only in the x and y 


coordinates. In this new parabolic form of the wave equation, the source term can be 


expressed as 


J(%, 0) =-ec> B,6° [% -F(0)]. {3.21} 


In this expression, /; (t) is the position of the ith electron at time ¢ and f, is the 


transverse velocity discussed in Equation {3.5}. For a helical undulator producing 


circularly polarized light, 7 = (—i,1,0), thus 


B, =Re [Sie {3.22} 
y 


Invoking the method of characteristics to introduce the coordinate u=z-—ct 
(which follows the light through the evolution), the derivatives in Equation {3.20} are 


transformed as follows: 
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Oz Ou {3.23} 
0 ome) 
Ot Ou Ot 
Recalling z as defined in Equation {3.1}, Equation {3.23} becomes 
. -ig 
[vi + ae E(%,t) = —4mieKk }o t), (3.243 
T ¥ 


where the notation <expression> indicates an average of the expression over all electrons 


and the electron particle beam density in a small volume element dV is 


p(%,t) = | ¥6 [FO] AV. {3.25} 


To simplify presentation (and mathematics), we utilize the definition of 
dimensionless laser field amplitude from Equation {3.14} and introduce the 


dimensionless current density 





__ 8Np(aeKL) 

jee {3.26} 
ymc 

Values of 7 <7 are “low gain” cases, while 7 >> is the “high gain” regime; as with 


the values of a, this FEL theory is valid for all values of. 


Transforming transverse coordinates to dimensionless form by x =x.J/k/(2L) 


and y= y,/k/(2L), the final form of the FEL wave equation is obtained: 


ix, 0 2 ay 


The modified Laplacian symbol Ve is used here to denote differentiation with respect to 
the normalized coordinates x and y. The terms of Equation {3.27} involving this 
spatial derivative describe the evolution of the optical field due to diffraction effects. 
Terms involving the time derivative describe its evolution due to gain from the electron 
beam. The term < je’ > is a measure of how strongly bunched the electrons are, 
showing that this portion of the dynamics have a significant effect on optical output. It is 


worth noting also that in Equation {3.27}, the evolution of a depends on electron 
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distribution ¢ , while in {3.14} the evolution of ¢ depends on |a|; thus, the FEL wave 
and pendulum equations form a feedback loop. Bunched electrons amplify light better, 
while the optical field’s coherent evolution promotes electron bunching. The gain from 
this process only grows for so long, however, as electrons will eventually overbunch, at 


which point the FEL reaches saturation. 


C: UNDULATOR TAPER 


In high-gain FELs, saturation can be reached before the end of the undulator, 
lowering overall FEL extraction (see below). This can be overcome by altering the 
resonance condition, Equation {3.11}, along the length of the undulator. While it is 
possible to do this by changing the undulator period 4, , it is more common to alter K. 
This typically involves changing the gap between undulator poles in order to change the 
magnetic field B, so the process is usually referred to as “tapering” the undulator. The 


effect is to add an additional phase acceleration term to the pendulum equation, giving 








2 
a > = 5 +|a|cos(¢ +), where 
dt 
ore {3.28} 
0 =-4nN ( } 
1+K°\ K 


The quantity 6 is referred to as the taper rate[4]. For taper rates that are too large 


(5>|a 





), the electrons are accelerated away from resonance before they are able to 


bunch, limiting energy transfer. Rates that are too low (6 < 4,Jha| ) will provide only a 


weak effect, and do not enhance the FEL’s natural extraction. Since oscillators typically 
have undulators too short to face this problem, tapering is much more commonly applied 


to FEL amplifiers[5]. 


D. EXTRACTION 


Up to this point, references to “efficiency” have applied to the overall system. To 
examine specifically the efficiency of the energy exchange between the electron beam 
and optical field, one may define the “extraction” as the relative amount of energy 


transferred from the electron beam to the optical field on a single pass, that is, 
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_ Ey, (AE) _ (Ar) 
ce eae leer a {3.29} 





where E,,, is the energy in the optical field, £, is that in the electron beam, and 


AE = Aymc’ is the change in the electron beam energy over the length of the undulator. 


Since the theoretical FEL gain curve has a width of Av = 27, it follows from Equation 


{3.13} that 


A Av 1 
Nmax a, ( 7) = Sie {3.30} 
y 4nN 2N 
Equation {3.30}, however, applies only to untapered undulators. Introducing taper 
changes the expression to [5] 
_jL 


3.31 
8aN j 


oO 


E. BEAM QUALITY 


The derivations of the FEL pendulum Equation, {3.14}, and the FEL wave 
Equation, {3.27}, assumed that the optical fields were plane waves of infinite transverse 
extent. This assumption remains valid mathematically so long as the electron beam 
remains close enough to the center of the optical field and diffraction is small. This 
places a condition both on the beam’s radius and on the angular spread of the electrons. 
Neither well-collimated beams of large radius nor highly divergent beams of small initial 
radius will satisfy the assumptions of Parts A and B. One appropriate measure of beam 


quality, then, is the transverse emittance, defined as 


£=1,0,, {3.32} 


where 7, is the root-mean-square (rms) radius of the electron beam and @, is its rms 


angular spread, as illustrated in Figure 4. 
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a) Well-Collimated Beam b) Highly Divergent Beam 





Figure 4. Beam Quality Dimensions 


An alternate measure, which has the advantage of being conserved as the beam is 


accelerated, is the normalized transverse emittance, defined as 


6 = ya {3.33} 


Equation {3.13} shows that changes in the phase velocity are proportional to the 
relative change in Lorentz factor, dy/y. This is nonzero for a relativistic beam from a 
particle accelerator because the electrons do not all have the same kinetic energy, and is 
referred to as the “energy spread.” It is another important measure of the beam’s quality, 
and it influences the FEL interaction, as well as the beam transport devices (e.g., focusing 


and defocusing elements). 
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IV. SIMULATION MODELS 


A. FEL3D/FEL4D 


Faculty and students at the Naval Postgraduate School have been engaged in 
theoretical research and numerical modeling of FELs for over twenty years. Their efforts 
have supported the FEL research community’s full spectrum of activity—from design 
support for new experiments to refining the understanding of experimental results. Over 
time, the Naval Postgraduate School (NPS) has developed a suite of simulations for both 
amplifiers and oscillators that have been benchmarked against both theoretical formulas 
and experimental results. The suite includes four types of programs: single-mode, 


longitudinal multimode, transverse multimode and four-dimensional [5]. 


The single-mode simulations follow the phase-space evolution of many electrons 
and one optical mode (an infinite plane wave) along the undulator, using the pendulum 
and wave Equations ({3.14} and {3.27}) to self-consistently determine the behavior of 
both the electrons and the light. They can be used to independently model certain 
physical effects, allowing for general understanding of the origins and consequences of 
those effects. They are small and computationally inexpensive, easy to understand, and 


provide a solid foundation from which to build more comprehensive simulations. 


The longitudinal multimode simulations are used to examine the evolution of 
longitudinal modes in the optical field. The optical amplitude and phase are allowed to 
vary as a function of z in each longitudinal “slice” of the optical field. Longitudinal 
modes play an important role in some FELs due to the phenomenon of “pulse slippage.” 
When the pulses are short, the difference in velocities between electron and optical pulses 


can significantly affect the manner in which the light is amplified. 


The transverse multimode simulations, collectively referred to as FEL3D, break 
from the assumption that the optical fields are plane waves. They use the pendulum and 
wave equations to self-consistently follow the evolution of the electron beam and the 
optical field in x, y, and t. They also allow for axial asymmetries, such as separate 
emittances in the two transverse planes, astigmatism in oscillator mirrors, and shifts or 


tilts of either the electron beam or the oscillator mirrors. These effects can also have a 
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significant impact on FEL gain and extraction. Both the transverse and longitudinal 
multimode programs require only moderate computing resources, and run on personal 


computers, though their efficiency is increased when multiple processors are available. 


To provide a more realistic simulation by following both longitudinal and 
transverse effects, NPS has recently developed the FEL4D programs. These model the 
electrons and the light in x, y, and z while allowing them to evolve in time. These 
simulations require significant processing time and large amounts of memory, and are 


used with the dedicated resources of a cluster computer. 


B. GINGER 


The history of the code Ginger began with the code Fred, a two-dimensional FEL 
amplifier simulation code developed at Lawrence Livermore National Laboratory 
(LLNL) some three decades ago. The early spinoffs from Fred included a harmonic 
simulation called Nutmeg and a three-dimensional (i.e., multimode) program dubbed 


Fred3D. 


Researchers at LLNL and Lawrence Berkeley National Laboratory (LBNL) 
created Ginger in the mid-1980s primarily to examine the physics and consequences of 
sideband growth in single-pass FEL amplifiers [6]. Ginger is polychromatic, time- 
dependent and multidimensional. It has grown over the years to include many optional 
features, such as a shot-noise startup model, configurations for complicated FEL types 
(e.g., harmonic cascade FELs), and higher-harmonic emission models. Historically, it 
has been limited to axisymmetric (7, z, tf) geometry for fields, but recent modifications 
include optional “slab” geometry (y, z, ¢, infinite x). The simulation can also still be run 


in a non-time-dependent (3D) configuration, referred to as “Fred mode.” 


Ginger and its postprocessor are both flexible and powerful. Both are written in 
Fortran90 and originally designed for use with Unix, but have been successfully ported 
to various platforms. The simulations explored for this thesis operate efficiently on a 
desktop computer; larger and more complicated problems can take advantage of the 


multiprocessor version of the program for use on cluster computers. 
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C. FEL3D CASE STUDY: BEAM QUALITY IN A GENERIC FEL 
OSCILLATOR 


Table 1 represents a parameter set for a generic FEL oscillator. These values do 
not correspond to any given machine, but are generally illustrative of realistic values and, 


thus, useful as a case with which to explore a simulation code’s operation. 


Energy: E,, = 120 MeV 
Radius: rp, =0.1 mm 
Bunch charge: q=1nC 
Bunch duration: th=2 ps 

Trans. emittance: & = 10 mm-mrad 
Long. emittance: & = 200 keV-ps 
Peak current: Teak = 500 A 
Energy Spread: Ay/y =0.2% 
Period: 

Length: 

Gap 

Number of periods: 

Undulator parameter: 


Wavelength: 
Length: 
Rayleigh length: 
Quality factor: 





Table 1. Generic FEL Oscillator Parameters 


Most of these quantities are self-explanatory or have been discussed before. 
Longitudinal emittance is an electron beam quality measure that is proportional to energy 
spread. The undulator gap g is the spacing between magnet poles, determining the 
magnetic field strength and, thus, the undulator parameter K. The Rayleigh length Zp is 
the distance from the optical beam waist over which the light will propagate before its 
cross-sectional area doubles. The optical wavelength value shown is the resonant 
wavelength of the FEL, determined from beam energy and undulator parameters by 


Equation {3.11} (recall that E, =ymc’). The quality factor Q, represents the losses per 


pass in the optical resonator, with low Q, values corresponding to high losses. In a 


conventional laser, this includes in-cavity losses such as heating in the lasing medium. 
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For an FEL oscillator, it is used to indicate the amount of power outcoupled through the 
resonator’s partially transmissive mirror. Ideally, it is related to the outcoupling mirror’s 
reflectivity R by 
1 
Pie {4.1} 
QO, 
The output of FEL3D is mostly graphical, though specific numeric quantities can 
be extracted with simple Unix commands. Figure 5 is part of the output for the case 


outlined in Table 1, depicting the +=Oslice of the optical field after 300 passes. The 


optical field amplitude a(0,¥,7) is shown, with lighter blue meaning more intensity. 


The thin cyan lines correspond to the e”' points relative to the peak amplitude. Spaced 
along the z axis are sample electrons (red dots)—the fact that they all fall within the inner 
portion of the optical field indicates that the beam and field were well-matched. Also 
depicted at c=0 and r=! are line-plots of the optical mode shape, with magenta 
corresponding to the cavity fundamental mode for reference and yellow to the mode 
observed in the simulation. Some deviation between the two is typical, especially at 


tT =1, and the “hollow” (depressed-peak) mode shown here is not uncommon. 





Figure 5. Generic Oscillator: Centerline Optical Field and Electron Beam (FEL3D) 





0 n 300 


Figure 6. Generic Ocillator: Power and Gain Evolution (FEL3D) 
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Another portion of the output is displayed in Figure 6. This depicts the evolution 
of outcoupled optical power and optical gain over 300 passes, in normalized units. Both 
quantities exhibit transient fluctuations for the first 100 or so passes, then reach steady- 


state. 


Several other figures are produced for a given case in FEL3D and can be 
reviewed either to check simulation viability or to look for particular behaviors. Often, 
however, the goal of a specific design (i.e., parameter set) is a particular quantity, usually 
the output power or extraction. To demonstrate the FEL3D code and gain a general 
insight into the behavior of FEL oscillators, we conducted a parameter study to determine 
the extraction against simultaneous variation of transverse emittance and energy spread. 


All other parameters used the values in Table 1. 
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Figure 7. Generic Oscillator: Beam Quality Parameter Study (FEL3D) 


The results of this study are shown in Figure 7. The results point to the general 
conclusion that increasing either emittance or energy spread reduces extraction. This is 
somewhat intuitive, but studies such as these can quantify the dependence. For these 
parameters, the peak extraction is approximately 3%, a strong but plausible value. The 


two red contours on this graph correspond to 2% and 1% extraction. Staying above half 
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the peak extraction requires simultaneously keeping the emittance under 30mm-mrad and 
the energy spread under 0.4%. Both requirements are reasonable values by current 


standards. 


D. GINGER CASE STUDY: JLAB’S 14KW EXPERIMENT 


In October 2006, Thomas Jefferson National Laboratory (JLab) broke their record 
for FEL average output power by achieving 14.2kW in continuous-wave oscillator 


operation. The parameters for that experiment are summarized in Table 2. 


Energy: E, = 150 MeV 
Radius: rp = 0.17 mm 
Bunch charge: q=114pC 
Bunch duration: t,=0.4 ps 
Peak current: Teak = 285 A 
Avg current: Lave = 8.5mMA 
Trans. emittance: &: = 8 mm-mrad 
Long. emittance: & = 80 keV-ps 
Energy Spread: Ay/y =0.2% 
Period: Ag=5.5 cm 
Length: L=1.65m 
Number of periods: N=30 
Undulator parameter: Kims = 1.36 
Length: 

Wavelength: 

Rayleigh length: 

Waist radius: 

Mode radius at mirrors: 

Quality factor: 





Table 2. JLab 14kW Parameters 


These parameters were used to build familiarity with Ginger just as the generic 
FEL oscillator parameters were used to build familiarity with FEL3D. A brief 
comparison shows that the numbers in Tables 1 and 2 are of similar magnitude. More 
importantly, Table 2 represents a known experiment with known results, which 


furthermore has been previously simulated by NPS. 
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The output of Ginger’s dedicated postprocessor is a series of graphs. Most of 
them track different quantities than the output of FEL3D, reflecting the difference in the 
programs’ origins. Gain and power evolution, however, are important to any FEL and 
appear in both codes’ output. For the case outlined in Table 2, the gain evolution is 


shown in Figure 8. 
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Figure 8. JLab 14kW: Gain Evolution (Ginger) 


As with the previous example, the simulation shows significant transients over the 
first one-third of the evolution, then settles into steady-state. The steady-state gain value 
appears to be approximately 25%. The gain per pass G required for steady-state 
operation is given by 

1 


reer {4.2} 


or 25% for Q, =5. This is a simple check that the simulation is behaving as expected 
and increases confidence in its results. 


The power evolution is depicted in Figure 9. The transients are minimal—a 
simple ramp-up—and the FEL reaches steady-state around pass 60, just as the gain did. 


This is another check on behavior meeting expectations. The values depicted on this 
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graph are of peak power at the exit of the undulator (P ). To convert to average 


outcoupled power (P,,,), one must know the pulse duration f, the pulse repetition 


frequency Q and quality factor Q,. The formula relating them is 


ie = Es (2) : {4.3} 
QO, 


For this experiment, the quality factor and pulse duration are given in Table 2, while the 


PRF is approximately 75MHz. Using 1.7GW for P_, yields P,, ~10.2kW. 


ut 
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Figure 9. JLab 14kW: Power Evolution (Ginger) 


While this result is approximately 27% below the known experimental results, it 
is of the correct order of magnitude. Furthermore, previously conducted NPS 3D 
simulations predicted nearly 20kW, a 40% overestimation. The conclusion drawn at the 
time was that pulse slippage effects were present in the experiment, which would reduce 


the predicted output power. 
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E. COMPARATIVE CASE STUDY: BEAM QUALITY IN A GENERIC FEL 
AMPLIFIER 


To compare the results of FEL3D and Ginger, we took the case of a generic FEL 
amplifier, using the parameters in Table 3. The values chosen match closely with the 
generic FEL oscillator parameters from Part C by design. The most notable difference is 
the significant increase in undulator length, typical for amplifier designs. Additionally, 
the optical cavity parameters have been replaced by seed laser specifications. The seed 
wavelength, which will also be the output wavelength, is determined entirely by the type 
of seed laser used and therefore typically driven by the needs of the application. In this 


case, the value was chosen to match the output wavelength of the generic oscillator. 


Energy: Ep =120 MeV 
Radius: rp =0.3 mm 
Bunch charge: q=l1 nC 

Bunch Duration: ty =2 ps 

Trans. Emit: €, = 10 mm-mrad 
Long. emit: €, = 200 keV-ps 
Peak Current: Tpeak = 500 A 
Energy Spread: Ay/y =0.2% 


Period: 

Length: 

Gap 

Number of periods: 
Undulator parameter: 


Wavelength: A=1.12 um 
Pulse Duration: top: = 4 ps 
Average power: Pon = 1 W 
Peak power: Preak = 330 kW 





Table 3. | Generic FEL Amplifier Parameters 


Using FEL3D, we conducted a parameter study on beam quality. As before, we 
varied transverse emittance and energy spread simultaneously and noted output 
extraction. For this version of the study, we included a taper rate of AB/B=5% over 


the last meter of the undulator. The results are shown in Figure 10. 
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Figure 10. Generic Amplifier: Beam Quality Parameter Study (FEL3D) 


Peak extraction for this case was somewhat lower than for the oscillator, near two 
percent. This is a normal and desirable value by current design standards. To stay above 
half this peak value, it was observed that the design must maintain transverse emittance 
below 10mm-mrad and energy spread below 0.25%. Recall that the corresponding limits 
for the oscillator were 30mm-mrad and 0.4%. This amplifier, then, is shown to be more 
sensitive to beam quality than a similar oscillator. This result is consistent with known 
experiments and the underlying physics. Put simply, amplifiers have significantly longer 
undulators; therefore, any degradation the beam quality introduces to the interaction is 
compounded over a longer time and has a greater effect. In an oscillator, each pass an 
optical pulse makes through the resonator is a new interaction with new electrons, so 


beam quality effects are essentially “reset” each time. 


For the code comparison, we examined the variation of emittance and energy 
spread separately. Figure 11 shows the results of varying emittance while holding energy 


spread constant at 0.2%. 
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Figure 11. Generic Amplifier: Extraction vs. Emittance (Comparative) 


The agreement between the codes is close. The only notable discrepancy occurs for the 
minimum emittance value. Neither simulation included taper in this case; the inset figure 
displays the power evolution from FEL3D, which shows clearly that the interaction 


reached saturation before the end of the undulator. 


Figure 12 shows the results of varying energy spread while holding emittance at 
10mm-mrad. Again, the codes agree very well until low values of energy spread. The 


inset indicates the same reasons as above for this discrepancy. 
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Figure 12. Generic Amplifier: Extraction vs. Energy Spread (Comparative) 
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F. CONCLUSIONS 


The two simulation codes considered—FEL3D and Ginger in its 3D mode— 
appear consistent. Some variation in the detailed results should be expected due to 
internal differences in the simulations. The strongest disparity seen in our examination 
came from the Jefferson Lab 14kW results, which Ginger underestimated and FEL3D 
overestimated. Both results, however, were of the correct order, and both discrepancies 
are likely due to the lack of pulse-slippage effects in the models. In direct comparison, 
the two codes agreed very well on beam-quality effects in a generic FEL amplifier. The 
discrepancies seen here show that the NPS models reached saturation early and therefore 


exhibited a lower extraction than they would have with a taper. 
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V. ELECTROSTATIC CATHODE TEST CELL DESIGN 


As discussed in Chapter II, design of an injector is important to overall design of a 
free electron laser (FEL). Characterization of the cathode used is important to the design 
of any injector. Normally, since injectors are generally not built to support more than one 
cathode type, this is performed on a “one-off” basis for the cathodes used by a given 
FEL. Additionally, replacing a cathode is often an involved task, undertaken only when 


the one in use has been spent or damaged. 


A tool the FEL community finds useful is the cathode test cell: a device that is 
not a full-fledged injector but can be used to test and characterize different cathodes 
under quasi-operational conditions. Due to the significant differences between the 
physical nature of and measurements required for thermionic, photoelectric, and field 
emission cathodes, the test cell design presented in this thesis focuses only on one type, 
namely field emission (FE). Existing FE test cells include the one used to explore 
diamond field-emitter array (DFEA) cathodes at Vanderbilt University (Figure 13). 
Typically, test cells achieve high fields with low voltages using small anode-cathode 
spacing (“gap”). At higher voltages, the small gap and gridded anode can result in 
subjecting the device and the cathode to significant damage from vacuum arcing and ion 
blowback. 


Collector: 
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Figure 13. Vanderbilt’s DFEA Test Stand [From 7] 
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The general goal for the design presented in this chapter is for a higher-voltage, 
larger-gap, non-gridded test cell. Building such a device will allow NPS researchers to 
explore the properties of FE cathodes at higher voltages and gain general experience with 
FE cathodes for possible use in their future FEL experimental facility. The design work 
discussed here presents preliminary considerations—building a foundation upon which 


future research will build, until a final design is selected and the system procured. 


A. FIELD EMISSION THEORY 


Cold (i.e., not thermionic) field electron emission can occur at the surface of a 
bulk material when strong static external electric fields are applied. It arises from the 
phenomenon of quantum tunneling and its explanation in the 1920s was one of the early 
triumphs in the then-emerging field of quantum mechanics. (It is also important due to 
its many practical applications, included devices such as the test cell under discussion.) 
The current density J emitted in the presence of applied field E is given by the Fowler- 


Nordheim equation, a simplified form of which is [8] 
J = AB’ /, {5.1} 


The constants A (in units of A/V”) and B (in units of V /m) are related to the material’s 


work function g and Fermi energy ¢, by 


a 6.2-10° le, 
pte, V@ {5.2} 
B=(68-10°)p” 


Several other, more complicated, forms of the Fowler-Nordheim equation exist, since 
several methods have been used to derive the fundamental form of the coefficients A and 
B. Equation {5.1} is presented in this form to illustrate the most salient points of the 
formula: inverse exponential dependence on £ and variation with material properties. 


The fundamental principle is that a larger field will yield more current. 


When the surface of the bulk material is not perfectly planar, any protrusions will 
cause local field enhancements (field lines will converge at their tips), effectively 
multiplying the applied field by a geometry-dependent factor at those points. This 


significantly lowers the overall field E required to reach a given value of current density 
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J. The field E as used in Equation {5.1} remains the pertinent value at the tips. This 
occurs naturally due to the microprotrusions (i.e., surface roughness) of bulk materials, 


but it can also be deliberately exploited by building in macroscopic protrusions. 


Diamond field emitter arrays (see Figure 14) take advantage of this technique. A 
thin layer of diamond is grown on a wafer of a suitable substrate such that the diamond 
forms into an array of pyramids. When placed in an external field, the pyramids serve as 
field-enhancement points and can emit electrons, which will be accelerated by the applied 
field and form beamlets. Arrays can be fairly dense (in terms of tips per unit area) to 
provide high current density, but this will cause space-charge forces to affect the 


beamlets’ propagation. 





Figure 14. Example of a DFEA [After 9] 


B. TEST CELL DESIGN REQUIREMENTS 


A cathode test cell must apply a given field to the cathode being tested and allows 
characterization or monitoring of the cathode’s performance. Proper characterization of 
an FE cathode involves measuring the current emitted as a function of field and time, and 
monitoring the beam shape, which will have an effect on beam transport and emittance. 
For DFEA cathodes, this is sometimes analyzed on a per-tip basis; for sparse enough 
arrays, the beam from each separate tip can be monitored. To support testing of multiple 
cathodes, the cathode holder should be easily removable. To facilitate varying the field 
strength FE, either the power supply voltage V must be variable, or the cathode assembly 
should be movable, allowing for adjustment of the anode-cathode gap d, or both. For the 
ideal case—where the anode and cathode are infinite parallel plates—the three quantities 


are related by 
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E=V/d. {5.3} 


Operation of a field-emission (FE) cathode requires a vacuum chamber—any gas 
present would be ionized by the electron beam and likely damage the cathode. Figure 
15 displays the insulating ceramic upon which this test cell design was based. Its metal 
flanges are separated by a ceramic break, allowing the components at either end to be 
electrically isolated from one another. The flanges are of a standard type (“conflat’), 
compatible with many common devices and fittings. Custom fittings can easily be 
ordered once the required configuration of measurement devices, power ports, vacuum 


system ports, etc., are fully determined. 


Le 





Figure 15. Ceramic Insulator with Conflat Flanges 


The electron beam’s size can be analyzed by placing a phosphor screen in the 
beam’s path and observing the spot(s) illuminated upon it. The phosphor screen can also 
serve as a charge collector, with current easily measured across a resistor placed between 


it and ground. The stability of the current measured can be monitored using a current 
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transformer placed around any isolated portion of the chamber through which the beam 
travels. Using the type of cylinder depicted in Figure 15, one can conceive of an 


arrangement such as shown in Figure 16. 


Here, the movable cathode assembly is fixed to an endcap that is bolted to the 
lower flange, allowing for easy removal. This end is held at a large negative voltage 
during operation. The annular anode is electrically coupled to the other endcap, which 
contains the phosphor screen (beam target) and is grounded during operation. Figure 16 


provides the basis for the preliminary test cell design presented in the following sections. 





Phosphor Screen 
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Figure 16. Electrostatic Cathode Test Cell Schematic 
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C; DESIGN TRADEOFFS AND PARAMETER SELECTION 


Experimentally, the emission predicted by Equation {5.1} is generally observable 
only once the applied field & has been raised to some minimum threshold value, 
dependent on the geometry and material. For the purposes of this thesis, the threshold 
field is considered to be 1OMV/m, and the maximum operating field is 20MV/m. These 
numbers are roughly characteristic of the DFEA cathodes to be studied at NPS. 


In addition to the insulating ceramic pictured above, NPS has on hand a fixed- 
voltage 100kV power supply suitable for use with this test stand. Field strength 
adjustment will therefore be made through adjustment of the anode-cathode gap. The 
pertinent tradeoff is between field strength and gap size: operating with a reduced the 
gap will achieve higher field but risk causing vacuum arcs, potentially damaging the test 


cell and cathode. 


This tradeoff also affects the placement of the phosphor screen. The simplest 
arrangement would be to fix the screen directly to the anode, but this creates significant 
risk of ion blowback as it is bombarded: ions excited off the screen would be positive 
and therefore be accelerated toward the cathode. This risk can be reduced by placing the 
screen past the anode in an approximately field-free region; any ions liberated from the 
screen will remain in its vicinity due to the lack of field instead of accelerating back to 
damage the cathode. This arrangement, however, requires an aperture in the anode, 
which will alter the field distribution and potentially defocus the beam. Some setups 
(including that in Figure 13) mitigate this by spanning the aperture with a grid or 
“pepperpot,” which mostly allows the electrons to pass, but keeps the field roughly 
uniform. This increases the risk of damaging the cathode and anode, however, albeit less 
so than with the direct-impingement arrangement. The grid’s focusing effects can also 
lead to emittance degradation. For this design, the aperture will remain ungridded; 
avoiding arcing and ion blowback is assumed to outweigh the effects of defocusing. This 


tradeoff remains a viable field of study for future design work. 


Though the electron beam will emerge from the cathode approximately collimated 
and perpendicular to the cathode surface, space-charge forces will tend to defocus the 


beam as it propagates. Techniques exist to counteract this, including the famous “Pierce 
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geometry” [10]. In a Pierce configuration, the metal surface on which the cathode rests is 
angled in the direction of beam travel, forcing the local field lines to point inward and 
help focus the beam. This is shown in Figure 17, where purple lines represent 
equipotentials and red arrows are local field vectors. The Pierce effect can help mitigate 
space-charge effects, but since it inherently enhances the field in certain regions, it 
increases the risk of breakdown. Pierce geometries also reduce the field at the cathode 
surface however (as is apparent in Figure 17), requiring a smaller gap for the same 
voltage. In this design, Pierce geometry is included at the cathode support but not 


analyzed in detail. The topic remains open for optimization in future studies. 
Cathode at -20 kV 


Grounded anode ees 





Figure 17. An Example of Pierce Geometry [After 11] 


The design’s major parameters are summarized in Table 4. The voltage value is 
fixed by available equipment, and the field strengths represent design goals; the range of 
gap values is the subject of study below, with initial estimates based on Equation {5.3}. 


Parameter 
Operating Voltage 100kV 


Anode-Cathode Gap 5-2.5mm 
Operating Field 10-20MV/m 





Table 4. | Guiding Design Parameters 
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D. THE DESIGN TOOL: POISSON SUPERFISH 


First developed at Los Alamos National Laboratory in the late 1970s, Poisson 
Superfish is a suite of programs used to calculate fields (electrostatic, magnetostatic or 
RF electromagnetic) in two-dimensional geometries. The component programs used for 
this design work include Automesh, which generates a triangular mesh (of variable 
spacing) fitted to the boundaries of the different regions of the geometry; Poisson, which 
calculates static fields by successive over-relaxation; the WSFPlot display tool, used to 
generate the diagrams shown below; and SF7, a postprocessor used to extract specific 


information from the field map (e.g., field strength along a specified line). 


Figure 18 displays an example of WSFPlot’s rendering of the region of the 
schematic indicated. Input geometries are specified in Cartesian coordinates, mapping 


x—r and y—z for cylindrically symmetric problems such as this. This is a common 


approach, and a flag within the input file will cause output tables to be presented as (7,z) 
instead of (x,y). Output files consist of a table representing the field map. WSF plot is 
used to overlay this on the problem geometry, and includes several useful display 


options, such as field arrows, equipotential contours and the mesh used. 
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Figure 18. Superfish Rendering of Test Cell Interior 
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The Superfish suite is a powerful tool for a wide variety of applications, and 
includes several additional programs not applicable to this design effort. There are also 
multiple postprocessors, including tabular display tools. For this thesis, however, all 


plots were generated in Matlab using data from SF’7. 


E. ANODE SHAPE 


The baseline concept for an anode in this design is a simple ring. A torus of metal 
with rectangular cross-section would be easily manufactured and serve the essential 
purposes of an anode. Figure 19 shows an initial design concept using this type of anode 
as a Close-up of the anode-cathode region. It is a portion of a cross-sectional drawing 
such as displayed in Figure 18 above; the left border is the test cell’s centerline. The red 
arrows depict the magnitude and direction of the field at given points, while the pink lines 


are equipotential contours. 





Figure 19. Simple Ring Anode 


The anode pictured is rectangular in cross-section, with its lower corner is filleted 


to avoid creating a field enhancement. This principle is also applied to the outer edge of 
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the cathode support stalk. As the field map shows, the field is elevated in that region 


even with the fillet; peak field value there should be monitored in future design iterations. 


The cathode is pictured in black, resting inside a divot in order to have Pierce 
focusing. The divot’s corners are not filleted, but do not appear to greatly perturb the 
field distribution. As expected, the field begins to defocus near the anode aperture, as 


evidenced by the sloping equipotentials. 


Figure 20 shows the second general design configuration considered, featuring an 
angled anode. Angling the lower surface should help alleviate the field concentration 
near the corner of the cathode support stalk. Angling the top surface helps mitigate 
damage to the anode by progressively increasing the distance between it and the electron 
beam. If the anode’s angle is greater than that of the diverging electron beam, no 
electrons will strike (and damage) the anode. Both of these features will affect field 
distribution; it is pertinent to ask whether they measurably affect field strength at the 
cathode. 





Figure 20. Angled Anode 
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F. GAP AND APERTURE SPACING 


Given a fixed voltage for the power supply, we studied two ways of altering the 
field distribution at the cathode: anode-cathode gap and anode aperture. The cathode’s 
own geometry has an effect as well, but was not studied as a part of this work. As a 
baseline, the initial configuration considered for each anode shape uses a lcm gap, as 
measured from the top of the cathode support stalk to the lowest point on the anode, and a 
lcm aperture. With a potential difference of 100kV, this gap value corresponds to a 
10MV/m gradient for infinite parallel plates (see Equation {5.3}). For the actual 
geometry, the results are significantly lower, as shown in Figure 21, due to Pierce effects 
and the anode aperture. Field magnitude values are shown as a function of radius out to 
5mm, the radius of the cathode divot. The cathode itself will have a radius of 
approximately 4mm. This baseline configuration data exhibits not only a failure to reach 
the minimum operating field but also a strong radial dependence in field strength. The 
fields are reasonably uniform only for the first 1-2mm, corresponding to 25% of cathode 
area. At a radius of 2mm, the axial field is reduced by 2.5% from its centerline value for 


both anode shapes. 


Gap=1cm, Aperture=1cm 
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Figure 21. Baseline Design 


The next case examined keeps the same gap value but narrows the aperture by 


half. If this showed improvement in field values, it could present significant advantages, 


4] 


chiefly by flattening out the field in the aperture. As Figure 22 shows, however, the 
improvement is slight, showing a peak field of 6.6MV/m vs. 6.4MV/m, and the strong 
radial dependence remains. As before, the simple ring anode provides somewhat more 
favorable fields, but in this case the disparity is a bit greater. Neither design appears 
overly promising, and the narrowed aperture introduces the same risks as a pepperpot 
design. A narrow aperture can even block a portion of the beam. While this is 
sometimes used in injectors to reduce beam emittance, it is not useful for measuring the 


cathode’s characteristics, and could damage the anode. 
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Figure 22. Narrow-aperture Design 


The next design iteration returns to the lcm aperture and halves the gap. The 
results, seen in Figure 23, show significant improvement. The peak field in this case 
actually exceeds the target value of 1OMV/m. The radial dependence of the field is 
somewhat lower but still significant, reaching the 2.5% reduction level around 2.5mm. 
The disparity between the two anode shapes has also narrowed, indicating that the angled 
anode may be a viable option at this gap spacing. While this case meets the minimum 


requirements, it is worth investigating yet smaller gap sizes. 


Decreasing the gap to 2.5mm produces the results shown in Figure 24. For this 
case, peak fields fall well within the operating range, and the angled anode performs 


almost as well as the simple ring. That result is favorable since such a narrow gap carries 
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Figure 23. Intermediate-gap Design 
increased risk of arcing and the angled anode may be vital to lowering fields outside the 


cathode region. Even more important a result is the field’s uniformity. The field strength 


only begins to drop off after about 3.5mm, meaning the field is uniform within 2.5% over 


75% of the cathode’s area. 


Gap=2.5mm, Aperture=1¢cm 
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Figure 24. Narrow-gap Design 
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G. RECOMMENDATIONS 


The conflat-cylinder insulator appears entirely feasible for testing cathodes. From 
a field-strength perspective, a 5mm gap is recommended as the maximum operating 
setting, and 2.5mm as the nominal setting, using a lcm or wider aperture. Future design 
work should examine the possibility of breakdown in the gap and whether an angled 
anode — perfectly suitable for achieving the required fields at the cathode—would provide 


mitigation. 
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